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I.  INTRODUCTION 


Although  accelerated  temperature  stressing  of  semiconductor  devices  is  an 

accepted  and  widely  employed  technique  for  assessing  the  reliability  of 

i  2 

semiconductor  devices,  *  it  appears  that  a  statistical  error  analysis  of  this 
procedure  is  not  available.  The  purpose  of  this  work  is  to  partially  remedy 
this  deficiency.  The  goal  of  this  program  is  easily  stated:  A  number  of 
devices  taken  from  a  distribution  are  operated  at  several  elevated 
temperatures  to  induce  failure  in  all  devices  within  a  reasonable  time. 
Assuming  general  characteristics  of  the  device  failure  probability  density 
function  (pdf)  and  its  temperature-dependence,  we  estimate  the  expected 
cumulative  failure  function  (cff)  for  devices  in  normal  operation.  By 
estimated  cff,  we  mean  our  best  estimate,  based  on  statistical  inference,  of 
the  average  probability  of  a  random  device  (taken  from  the  same  distribution 
but  operated  at  a  normal  temperature)  failing  as  a  function  of  time. 

Section  II  contains  a  brief  review  of  the  general  mathematical  formalism 
usually  employed  in  semiconductor  reliability  discussions.^  Three  failure 
pdf's  ot  particular  usefulness  to  this  analysis — exponential,  normal,  and 
lognormal — are  discussed.  Our  particular  interest  is  in  highly  complex 
electronic  systems  intended  for  long-life  space  applications.  In  this 
application  and  for  similar  ground-based  systems,  it  is  not  the  expected 
failure  probability  density,  at  times  comparable  to  the  median  lifetime  ("10^ 

Q 

to  10  hr)  of  the  device  that  is  important,  but  the  cff,  at  times  orders  of 
magnitude  less,  at  times  comparable  to  the  desired  system  useful  life  (“10^  to 
10^  hr).  A  simplified  cost  analysis  to  justify  this  assertion  is  included  at 
the  end  of  Section  II.  A  brief  review  of  accelerated  temperature  aging  is 
presented,  and  the  assumptions  concerning  the  general  characteristics  of  the 
failure  pdf,  which  are  fundamental  to  this  analysis,  are  emphasized. 

This  analysis  is  carried  out  in  several  steps,  each  more  difficult  than 
the  preceding.  We  first  consider  the  case  of  operating  the  devices  to  failure 
at  a  single  temperature.  On  the  basis  of  the  experimental  observed  failure 
times  and  the  assumption  of  a  lognormal  pdf,  we  estimate  the  lifetime 
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distribution  of  the  original  population  from  the  failed  sample  distribution. 
The  maximum-likelihood  method,  combined  with  the  chi-squares  and  the  Student-t 
distributions.  Is  employed  to  estimate  the  median  lifetime,  logarithmic 
variance,  and  their  associated  distributions,  as  well  as  the  average 
cumulative  failure  function. 

This  procedure  Is  repeated  for  a  generalized  accelerated-temperature¬ 
stressing  experiment,  and  these  results  are  first  used  to  analyze  a  simple 
two-temperature  stress  test  because  the  physical  understanding  is  not  obscured 
by  the  mathematics  and  also  because  such  a  program  is  statistically  very 
efficient  within  the  constraints  of  the  accelerated  temperature  range  or 
equivalently  within  the  constraint  of  limited  test  time.  Various  other  sample 
distributions  within  the  elevated  temperature  range  are  considered,  and  It  Is 
shown  that  samples  near  the  center  of  the  accelerated-temperature  region 
contribute  (statistically)  inefficiently  to  the  accuracy  of  the  extrapolated 
estimates.  The  final  result  of  this  analysis  Is  an  estimate  of  the  average 
cumulative  failure  function  of  a  device  operated  at  the  actual  operating 
temperature.  This  estimate  is  based  on  the  sample  failure  distribution  of  a 
thermally  accelerated  aging  test  with  no  prior  knowledge  of  either  the  median 
lifetime  or  its  logarithmic  variance  assumed. 


II.  BACKGROUND 


The  reliability  of  a  given  device  is  conveniently  described  by  its 
failure  pdf,  which,  by  convention  is  designated  f(x),  where  x  is  a  random 
variable.  The  differential  probability  that  a  randomly  chosen  device  will 
fail  between  x  and  x  +  Ax  is  f(x)Ax.  The  cff,  F(t),  is  the  probability  that  a 
given  device  will  have  failed  before  t  and  is. 


F(t) 


f (x)dx 


(1) 


The  survival  probability  is  one  minus  the  cff,  when  one  assumes,  as  is  done  in 
the  present  analysis,  that  a  device  can  only  be  in  one  of  two  states,  i.e., 
fully  operational  or  failed.  The  failure  rate,  X(t),  is 


X(t) 


f(t) 

[I  -  F(t) ] 


(2) 


A  given  pdf  can  be  characterized  by  many  different  parameters,  but  only 
three  of  these  will  be  employed  in  this  analysis:  mean  lifetime,  median 
lifetime,  and  the  variance  of  the  distribution.  The  mean  lifetime  for  a  given 
pdf  is  designated 


t  - 


tf(t)dt 


(3) 


i 
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and  represents  approximately  the  average  life  expectancy  of  a  large  number  of 
devices  taken  for  the  distribution  under  discussion.  The  variance  is  a 
measure  of  the  "spread"  of  the  lifetimes  and  is  defined  as 


■  (  (t  -  t)^f(t)dt 

*'o 


(4) 


In  this  case,  the  variance  is  defined  with  respect  to  "linear"  time  as  the 
random  variable.  Although  this  is  the  usual  designation,  the  choice  of  the 
random  variable  is  not  unique.  With  semiconductor  devices,  it  is  convenient 
to  specify  the  variance  with  respect  to  the  logarithm  of  operating  time  as  the 
random  variable.  The  median  lifetime  is  similar  to  the  mean  lifetime  except 
that  it  represents  the  time  at  which  the  cff  is  equal  to  0.5.  If  a  large 
sample  of  devices  were  operated  to  failure,  the  median  lifetime  or  log 
lifetime  is  approximately  equal  to  the  time  at  which  one-half  of  the  devices 
have  failed. 

The  general  properties  of  three  failure  pdf  which  are  of  particular 
interest  to  the  present  analysis  are  briefly  reviewed. 

The  exponential  failure  pdf  is 


f(t)  *  X  exp(-Xt) 


(5) 


where  X  is  a  constant.  For  any  true  pdf,  we  have 


f  f(t)dt  -  1 


(6) 
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The  cff  assumes  a  simple  form 


F(t)  =  1  -  exp(-Xt) 


(7) 


”  Xt  for  Xt  «  1 


(8) 


The  failure  rate  is  constant  and  equal  to  X.  The  mean  lifetime  and  its 
standard  deviation  are  both  equal  to  the  reciprocal  of  the  constant  failure 
rate. 

The  normal  failure  pdf  is 


f(t) 


o(2ir) 
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exp 


(9) 


9 

where  g  is  the  mean  lifetime  and  a  the  variance.  Unlike  the  exponential 
distribution,  the  median  lifetime  is  equal  to  the  mean  lifetime.  The  cff 


F(t) 


/ 


o(2tt) 


1/2 


exp[~ 


(10) 


cannot  be  expressed  in  a  closed  form  and  must  be  evaluated  by  using  a 
tabulation  of  the  standard  normal  distribution 


9 


F(a) 


(2it) 
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/(a-p)/a  .  « 

e*p(-  \  y  )dy 


(ID 


Probability  graph  paper  (e.g.,  K  &  E/46  8003)  is  particularly  useful  when 
working  with  normal  distributions.  When  the  cff  is  plotted  on  this  graph 
paper,  a  straight  line  results.  Both  the  median  and  the  mean  lifetime  occur 
at  cff  *  0.5  (i.e.,  50%)  and 


a  «  p  - 


16 


(12) 


where  tjg  is  the  time  corresponding  to  F(t)  -  0.16.  The  normal  pdf  is 
symmetric  about  p,  and  approximately  68%  of  the  failures  occur  within  o  of  the 
median  lifetime. 

Under  carefully  controlled  conditions  of  elevated  temperature  and  applied 
bias,  the  failure  pdf  of  semiconductor  devices  is  found^’^  to  follow  a 
lognormal  pdf  expressed  as 


f(t) 


1 


to(2ir) 
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ex 


.  ln(t/xM)  2 

*[-  ^  n  ] 


(13) 


In  this  expression,  is  the  median  lifetime,  and  a  is  the  logarithmic 
standard  deviation.  The  lognormal  distribution  is  similar  to  the  normal  pdf 
except  that  t  ♦  In  t.  In  fact,  if  we  define 


x  ■  In  t 


(14) 


10 


we  have 


f(x) 


_ 1_ 

o(2x) 


172  exp[_ 


1/-X  -  u 
2 a 


(15) 


a  normal  distribution  employing  not  a  linear  time  scale  but  a  logarithmic  time 
scale.  By  working  In  logarithmic  time,  all  the  properties  of  the  normal 
distribution  apply  to  the  lognormal  pdf.  The  time  interval  [t,  t  +  At] 
becomes  [x,  x  +  Ax]  *  (In  t,  In  t  +  At/t].  The  cff  is 


F(x)  - 


/ 


1 


o(  2n) 
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ex 


p[-  -k2-^*1)2] 


(16) 


which,  when  plotted  on  probability  graph  paper  using  the  logarithm  of  time, 
results  in  a  straight  line.  This  straight  line  crosses  the  50%  probability  at 
x  *  u;  the  median  of  the  log  lifetime  distribution.  The  median  lifetime  Is, 
of  course. 


tm  *  exp  p 


(17) 


In  a  manner  similar  to  a  normal  distribution,  we  have  approximately 


°  "  M  “  X16 


(18) 
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where  Xjg  and  t|g  correspond  to  the  value  of  x  and  of  t  where  16%  of  the 
cumulative  failure  occurs.  Approximately  68%  of  the  failure  occurs  between 
xj£  and  xg^  (i.e.,  2o,  centered  around  the  median  lifetime). 

Expressions  for  the  real-time  (i.e.,  linear)  mean  lifetime  and  variance 
[i.e.,  Eqs.  (3)  and  (4)]  are  available  but  serve  no  useful  purpose  in  the 
present  analysis.  In  fact,  we  emphasize  that  knowledge  of  the  failure  pdf 
near  the  mean  lifetime  is  of  little  value  in  analyzing  the  reliability  of 
complex  systems  intended  for  long-life  space  applications  or  other  similarly 
complex  ground-based  systems.  Of  prime  Importance  is  the  failure  rate  at 
times  much  less  than  the  average  lifetime.  In  general,  with  a  system 
consisting  of  n  components,  one  desires  to  know  the  cff  near  the  100/n 
percentage  level.  For  example,  a  system  consisting  of  10^  identical 
components  in  series  would  be  expected  to  have  a  high  probability  of  failure 
at  a  time  when  the  average  device  cff  is  10~^. 

If  we  look  at  a  total  system  from  the  viewpoint  of  cost-effectiveness,  it 
is  not  the  failure  pdf  or  cff  that  one  is  interested  in,  but  rather  the 
integral  of  the  cff  over  the  service  life.  Consider  a  system  with  a  specified 
useful  life  of  t8«  Let  the  loss  in  some  monetary  units  for  the  systems  not 
working  in  the  time  interval  [t,  t  +  At]  be  R(t)At.  The  expected  loss  from 
failure  of  the  it^'  component  is 


AL  - 
s 


J  R(t)F1(t)dt 

o 


(19) 


where  F^(t)  is  the  cff  for  the  i^  component.  The  loss  rate,  of  course, 
depends  on  the  system-intended  employment  and  reliability.  Assuming  a 
constant  "usefulness"  throughout  the  system  lifetime,  the  loss  rate  is 
expected  to  be  proportional  to  the  survival  probability  of  the  system 
excluding  the  i^  component.  Assuming  a  high  survival  probability  for  the 
total  system,  the  loss  due  to  the  failure  of  any  component  is  independent  of 
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the  failure  of  other  components.  If  the  rate  of  loss  R(t)  can  be  reasonably 
represented  as  a  constant,  R,  we  estimate 


(20) 


or  the  total  estimated  loss  from  all  components,  assuming  EF^(t)  «  I  for  t  < 
t8,  is  the  summation  of  the  individual  losses  and  is 


L.  /%1(t)dt 


(21) 


For  an  exponential  pdf,  this  takes  a  particularly  simple  form  assuming 
A^t  «<  1,  and  we  write 

\  ■  T  *  h  '22> 

i.e.,  one  simply  sums  the  various  failure  rates.  For  a  lognormal 
distribution,  a  much  more  complex  situation  results  because  a  single  lognormal 
cff,  let  alone  the  sum  of  lognormal  cff's,  cannot  be  Integrated  to  yield  a 
closed  form  solution. 

Thus  far,  we  have  been  speaking  as  if  we  knew  the  failure  pdf 
accurately.  In  practice,  we  do  not  know  the  failure  pdf  exactly;  in  fact,  in 
many  cases,  we  do  not  even  know  the  functional  form  of  the  pdf.  Only  by 
sampling  from  the  distribution  can  we  estimate  the  correctness  of  an  assumed 
functional  form,  and  only  when  this  is  established  are  we  in  a  position  to 
estimate  the  parameters  associated  with  the  failure  pdf.  Knowing  the 
functional  form  of  the  failure  pdf,  one  estimates  the  expected  cff  from  the 
sample  distribution  and  not  the  exact  cff.  In  the  present  analysis,  it  is 
important  to  distinguish  between  the  parameters  that  characterize  an  exact  pdf 
and  the  corresponding  parameters  estimated  from  sampling. 
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III.  ACCELERATED  TEMPERATURE  TESTING 


Because  semiconductor  components  have  excellent  reliability,  it  is  very 
costly  to  adequately  determine  the  failure  pdf  under  actual  operating 
conditions.  This  excessive  cost  is  associated  with  the  large  sample  of 
devices  that  must  be  operated  for  long  test  times  to  adequately  predict  the 
failure  pdf  with  a  reasonable  confidence  level.  Even  if  the  cost  of  such  a 
reliability  test  is  not  a  constraint,  the  required  total  test  time  is  usually 
a  constraint.  If  one  desires  to  predict  the  reliability  of  an  electronic 
system  with  a  desired  system  life  of  tg,  one  must  determine  the  failure  rate 
at  times  comparable  to  tg,  which  for  space  applications  can  be  ten  years  or 
longer.  If  the  failure  pdf  is  exponential,  one  can  decrease  the  total  test 
time  by  increasing  the  sample  size,  since  in  this  special  case,  it  is  the 
total  operating  time  of  the  device  that  establishes  the  accuracy  of  predicting 
the  failure  pdf.  When  the  functional  form  of  the  failure  pdf  is  unknown,  it 
is  dangerous  to  assume  that  one  can  trade  sample  size  against  test  duration. 

Because  of  this  consideration,  one  often  resorts  to  an  accelerated 

temperature-stress  program  to  estimate  the  reliability  of  semiconductor 

components.  The  success  of  such  a  program  is  based  on  the  remarkable 

1  2 

experimental  observations  ’  that  for  a  single  failure  mechanism 

1.  The  failure  pdf  at  a  constant  temperature  and  applied  electrical 
stress  is  lognormal. 

2.  The  logarithmic  variance  is  independent  of  temperature. 

3.  The  median  failure  time  follows  an  Arrhenius  dependence  expressed  as 


To  exp(kf) 


where  xQ  and  AE,  the  activation  energy,  depend  on  the  electrical  stress  but 
not  on  temperature.  A  typical  accelerated  test  is  schematically  shown  in 
Fig.  1,  in  which  a  small  sample  of  devices  are  operated  to  complete  failure  at 
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Fig.  1.  Arrhenius  Presentation  of  Accelerated  Temperature  Aging.  This 

presentation  was  constructed  by  picking  10  random  lifetimes  at  three 
different  temperatures  from  a  lognormal  distribution  whose  parameters 
are  u  -  In  [10“®  exp  (0.794/kT)]  and  o  *  1.00.  The  exact  median 
lifetimes  at  the  various  temperatures  are  listed,  as  well  as  the 
sample  mean  lifetime  and  the  logarithmic  standard  deviation.  The 
solid  line  represents  the  Arrhenius  dependence  of  the  exact  median 
lifetime,  which  extrapolates  to  ■  10®  hr  at  an  operating 
temperature  of  60°C. 
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three  different  temperatures.  The  logarithms  (to  the  base  10)  of  the 
experimentally  measured  failure  time  are  plotted  as  a  function  of  the 
reciprocal  of  the  absolute  "junction"  temperature.  The  scatter  in  the  failure 
times  at  constant  junction  temperature  is  the  result  of  the  random 
distribution  of  failure  times  as  expressed  by  their  failure  pdf.  The  average 
increase  in  the  failure  time  with  increasing  reciprocal  absolute  junction 
temperature  is  because  their  median  lifetime  follows  an  Arrhenius  dependence. 

If  one  could  take  a  given  device  to  failure  at  an  accelerated 
temperature,  restore  the  device  to  its  original  condition,  and  then  repeat  the 
failure  at  several  different  temperatures,  the  log  of  the  observed  failure 
times  plotted  as  a  function  of  the  reciprocal  junction  temperature  would  be 
expected  to  follow  a  straight  line  characterized  by  a  specific  tq  and  AE. 
Repeating  this  experiment  with  other  devices  from  the  same  uniform  pdf  would 
result  in  a  series  of  parallel  lines,  indicating  constant  failure  activation 
energy.  The  vertical  shift  between  these  Arrhenius  dependences  is  caused  by 
the  scatter  in  the  pre-exponential  factor,  tq.  Since  the  pre-exponential 
factor  in  an  Arrhenius  dependence  is  relatively  temperature  insensitive,  the 
associated  logarithmic  variance  is  also  temperature  insensitive. 

Thermally  accelerated  aging  results  are  typically  analyzed  as  follows: 

The  failure  times  observed  at  a  given  temperature  are  plotted  on  probability 
paper  to  determine  whether  the  sample  can  reasonably  be  assumed  to  have  come 
from  a  distribution  whose  pdf  is  lognormal.  One's  confidence  in  being  able  to 
make  such  a  judgment,  of  course,  increases  with  the  sample  size  failed  at  any 
given  temperature.  If  the  failures  at  each  temperature  follow  a  lognormal 
pdf,  the  logarithmic  variance  at  each  accelerated  temperature  is  estimated  to 
determine  whether  it  is  temperature  independent.  As  before,  this  decision  is 
greatly  aided  by  a  large  number  of  failures  at  each  temperature.  The 
logarithmic  variance  usually  ranges  from  0.5  to  2.0.*  A  variance  of  less  than 


^Reasonable  care  must  be  exercised  in  this  type  of  analysis  because  of  rapid 
interchange  that  occurs  between  naperian  base  (natural)  logarithm  and 
logarithms  to  the  base  20.  The  natural  logarithm  is  used  for  o*  in  this  work. 


17 


1.0  is  characteristic  of  a  reasonably  mature  and  well-controlled  process, 
whereas  a  variance  of  greater  than  2.0  is  Indicative  of  an  immature  process 
and  is  normally  rejected  out  of  hand.*  Having  established  a  reasonably  low 
variance  that  is  supposedly  temperature  insensitive,  we  plot  the  results  as 
indicated  in  Fig.  1  to  decide  whether  the  failures  can  be  fitted  with  a  single 
activation  energy.  One's  confidence  in  making  such  a  judgment  once  again 
increases  with  sample  size.  If  it  is  decided  that  a  single  activation  energy 
is  justified,  constants  and  AE  are  found  that  best  describe  the  data.  With 
these  constants,  and  the  logarithmic  variance,  one  can  then  determine  the  most 
appropriate  pdf  for  any  desired  operating  temperature. 

This  common  procedure  has  a  serious  flaw  that  should  be  emphasized.  The 
failure  pdf  is  estimated  as 


,  ,  In  t  -  2 

f (In  t)  -  - - J72  exp[-  - = - )  ] 

o( 2v)l,Z  i  o 


(24) 


—  *2 

where  Ujj  and  o  are  the  experimenter's  best  estimate  of  the  extrapolated 
median  log  lifetime  and  the  variance.  As  such,  it  estimates  the  "median”  cff 
but  not  the  "average"  cff,  which  can  be  very  different,  especially  for  times 
much  shorter  than  the  extrapolated  median  lifetime.  The  important  distinction 
between  "median"  and  "average"  estimated  cff  is  the  primary  motivation  of  this 
analysis  and  is  discussed  in  great  detail  later. 

Our  analysis  follows  the  standard  treatment  given  in  many  books  on  the 
theory  of  statistics.  The  authors'  preference  is  the  early  edition  of  Mood 
[3].  An  excellent  bibliography  on  the  theory  of  statistics  is  presented  in 
Ref.  [4].  An  extensive  discussion  of  the  lognormal  distribution  has  been 
provided  by  Altchison  and  Brown  15], 

Jordan  {61  presented  a  comprehensive  review  of  the  lognormal  failure 
distribution  to  the  analysis  of  semiconductor  component  failure.  His  analysis 
differed  from  the  present  work  in  that  he  confined  his  attention  to  estimating 
the  statistical  error  associated  with  the  median  lifetime  and  the  logarithmic 
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variance  from  a  sample  of  devices  all  aged  to  failure  at  a  constant  stress 
temperature.  Our  primary  intention  is  to  estimate  the  error  associated  with 
extrapolating  thermally  accelerated  failures  to  a  lower  operating  temperature. 

An  extensive  body  of  literature  on  the  statistical  analysis  of 
accelerated  stress  testing  has  been  contributed  by  the  General  Electric  (GE) 
Research  and  Development  Center  staff.  This  work  does  not  appear  to  be  widely 
discussed  within  the  semiconductor  device  reliability  community  as  is  evident 
in  that  it  is  not  referenced  in  {1],  [2],  and  [6].  Perhaps  the  explanation 
for  this  unfortunate  oversight  is  that  the  GE  group  did  not  employ  accelerated 
temperature  aging  of  semiconductor  device  results  as  a  vehicle  to  illustrate 
the  application  of  their  statistical  analysis.  Regardless  of  what  test 
results  were  used,  their  assumptions  (i.e.,  lognormal  pdf /Arrhenius-dependent 
median  log  lifetiwO  are  the  same  as  commonly  employed  in  semiconductor  device 
reliability  discuss ’ips,  and  their  work  is,  therefore,  directly  applicable. 

The  GE  analysis  is  conveniently  separated  into  two  parts;  complete  and 
censored  data.  Ir,  Nelson's  early  work  [7],  the  analysis  of  accelerated 
temperature  uglng  to  complete  failure  was  considered  in  greater  detail  than 
given  in  Ref.  [8].  His  main  interest  was  to  estimate  the  confidence  limits  at 
the  extrapolated  median  lifetime  and  the  associated  standard  deviation. 
However,  our  main  interest  is  to  estimate  the  probability  that  a  device  will 
fail  at  a  time  much  less  than  the  extrapolated  median  lifetime,  a  subject 
given  only  cursory  consideration  by  Nelson.  This  difference  in  emphasis  was 
motivated  by  the  final  application  of  the  results.  We  are  interested  in 
highly  complex  systems  where  the  failure  rates  at  times  much  less  than  the 
median  lifetime  are  important ^whereas  Nelson's  interest  appears  to  have  been 
in  systems  with  few  components,  and,  therefore,  failure  rates  near  the  median 
lifetime  are  emphasized.  Regardless  of  this  difference  in  emphasis  between 
the  present  work  and  the  previous  accelerated-aging  analysis,  there  are  no 
basic  differences  or  conflict  because  they  both  follow  the  standard  linear 
regression  analysis. 

In  more  recent  work  from  the  GE  group,  the  analysis  of  censored  data  is 
considered  [9].  An  example  of  censored  data  in  connection  with  accelerated- 


temperature-stress  testing  would  be  where  only  the  first  five  failure  times  at 
each  temperature  would  be  available  in  the  experiment  discussed  in  connection 
with  Fig.  1.  The  analysis  of  censored  data  is  much  more  complex  than  the 
problem  of  complete  failure  data.  The  subject  of  censored  data  is  usually 
described  in  more  advanced  probability  and  statistical  books,  such  as  Kendall 
and  Stuart's  second  volume  [10].  However,  the  Important  subject  of  the 
analysis  of  censored,  accelerated-aging  data  is  outside  our  present  objective. 

The  correctness  of  an  accelerated-temperature-aging  program  depends  on 
systematic  and  statistical  errors.  Systematic  errors  can  seriously  compromise 
the  effectiveness  of  an  accelerated  stress  program.  A  lower  activation  energy 
failure  mechanism  invalidates  extrapolating  the  higher  temperature  failures  to 
lower  temperatures.  Although  the  authors  recognized  that  a  thermally 
accelerated  stress  program  is  very  prone  to  systematic  errors,  they  have 
confined  this  analysis  to  statistical  errors.  In  this  context,  the  results  of 
this  analysis  will  provide  necessary  but,  of  course,  not  sufficient  criteria 
to  help  determine  whether  devices  from  a  particular  lot,  a  sample  of  which  was 
subjected  to  thermally  accelerated  aging,  should  be  installed  in  a  particular 
system.  The  conclusions  of  this  work  are  valid  only  under  the  following 
assvmptions: 

1.  The  failure  pdf  is  described  by  single  lognormal  distribution. 

2.  The  lognormal  variance  is  independent  of  temperature. 

3.  The  median  lifetime  follows  an  Arrhenius  dependence  whose  pre- 
exponential  factor  and  activation  energy  are  independent  of 
temperature. 

It  is  recognized  that  these  assumptions  apply  only  to  an  unbiased,  i.e., 
electrically,  as  opposed  to  statistically,  accelerated  stress  program  and  that 
the  lifetimes  should  be  determined  from  a  biased  accelerated-temperature 
stress  program.  The  above  assumptions  are  too  restrictive  in  the  sense  that 
if  the  bias  dependence  of  o,  xQ,  and  AE  are  kept  constant  during  the 
accelerated  stressing  and  equal  to  the  values  at  the  operational  level,  the 
results  of  this  analysis  would  also  apply  to  a  biased  accelerated  stress 
progrsm.  The  authors  have  chosen  not  to  include  the  biased  accelerated  stress 


assumption  not  because  they  feel  it  is  invalid,  but  because  they  wish  to 
emphasize  the  importance  of  carefully  examining  the  implication  of  electrical 
bias  during  accelerated  stress  testing.  Also,  the  assumption  of  a  single 
lognormal  pdf  rules  out  the  possibility  of  the  existence  of  a  small  percentage 
of  early  failures. * 
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IV.  STATISTICAL  ANALYSIS  OF  ACCELERATED  TEMPERATURE  AGING 


In  a  single-temperature  test,  n  devices  are  operated  until  they  all  fail. 
Based  on  the  n  different  failure  times  experimentally  observed,  one  desires  to 
estimate  the  population  distribution  from  this  sample  distibution.  Establish¬ 
ing  the  validity  of  an  assumed  functional  form  for  the  failure  pdf  is  outside 
the  scope  of  this  analysis  and  we  proceed  directly  to  estimating  the  lognormal 
parameters  from  the  experimentally  observed  failure  times,  t  ^ ,  t^,  T3•••  This 
will  be  done  by  the  commonly  employed  maximum  likelihood  method.  The  results 
of  this  analysis  of  a  single-failure  temperature  closely  follows  the  analysis 
of  Jordan^  and  is  included  here  both  for  completeness  and  for  its  usefulness 
in  introducing  the  more  complex  analysis  to  follow. 

The  experimental  observed  failure  times  are  converted  to  the  correspond¬ 
ing  logarithmic  X^,  X2,  X3...,  where  each  of  these  samples  is  from  a  distri¬ 
bution  whose  median  value  and  variance  we  wish  to  estimate.  In  the  maximum- 
likelihood  method,  one  designates  the  deviation  of  the  sample's  values  from 
the  median  value  and  forms  the  summation 

&  -  If(Xt  -  y)  (25) 


or 


In  !£  -  -  nln  [o(2n  - ^(Xj^  -  y  )2  (26) 

la 

where  the  summation  is  performed  over  the  complete  sample  distribution  and  is 
referred  to  as  the  log-likelihood  function.  This  function  contains  all  the 
information  available,  i.e.,  the  assumed  pdf  functional  form  and  the 
experimental  results,  and,  therefore,  represents  a  reasonable  starting  point 
in  any  statistical  analysis.  In  the  method  of  maximum  likelihood,  it  is 
assumed  that  reasonable  estimators  for  y  and  a  are  the  values  y  and  0,  which 
maximize  the  maximum-likelihood  function  or,  more  conveniently,  the  log- 
.  likelihood  function. 


Equating  to  zero  the  derivatives  of  the  log-likelihood  function  with 
respect  to  both  p  and  a  provides  the  desired  point  estimators,  which  are 


(27) 


and 


a2  -  £  £(Xt  -  p)2  (28) 

_2 

The  estimator  for  the  median  value  is  unbiased;  o  is  negatively  biased  but 
consistent,  i.e.,  consistent  in  that  in  the  limit  of  large  sample  size 
a+0.  The  best  estimator  of  a  is,  of  course, 

O2  -  ^—Y  I(*i  -  V)2  (29) 

If  we  know  a  exactly,  we  can  use  Eq.  (27)  and  the  self-reproductive 
property  of  normal  populations  to  estimate  the  distribution  of  the  median 
value  as 

om(2it)  M 


where 


(31) 


and  is  commonly  referred  to  as  the  standard  deviation  of  the  mean.  This 
equation  expresses  the  simple  truth  that  the  error  in  the  estimated  median 
value  follows  a  normal  distribution  whose  variance  decreases  linearly  with 
sample  size.  If  a  large  number  of  workers  were  to  repeat  this  experiment, 
their  composite  results  for  jT  would  follow  the  above  expression  except  that 
their  median  IT  would  be  nearly  equal  to  the  exact  p  because  of  the  large  total 
sample  size  Involved.  The  variance  in  the  reported  values  of  p  from  all  these 
workers  would  be  nearly  equal  to  the  variance  expressed  in  Eq.  (31). 
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The  confidence  in  one's  ability  to  predict  o  from  the  sample  distribution 
can  be  estimated  by  employing  the  well-known  chi-square  pdf.  For  a  normal 
distribution,  the  random  variable 


2  IO^-m)2 


has  the  chi-square  distribution  with  n-1  degrees  of  freedom.  The  probability 

2 

that  a  randomly  chosen  value  of  x  is  within  some  interval  (a,b)  is 


P{a  <  I(X1  -  p)2  <  b}  -  fbf  (x2.!  )dx  2_! 

n 


P(2L.  <  o2  <  SL  )2 
v  b  a 


where  6  can  be  found  from  cumulative  chi-square  distribution  tables  once  the 
desired  confidence  interval  is  specified. 

Knowing  neither  p  nor  a,  one  can  determine  the  expected  uncertainty  in 

p  by  employing  the  Student-t  distribution.  Most  standard  f.exvs  on  statistical 

analysis  prove  that  if  x  is  a  normally  distributed  random  variable  with  mean 
2  2 

p  and  variance  o  ,  if  x  has  a  chi-square  pdf  with  k  degrees  of  freedom,  and 
2 

if  x  and  x  are  independently  distributed,  the  function 

t  _  (x  -  p)/o 


[x2/*]1'2 


Our  reproducing  the  complex  derivation  of  the  chi-square  pdf,  or  of  the 
Student-t  pdf  to  be  shortly  Introduced,  serves  Jittle  purpose  in  our  primary 
goal.  We  desire  only  to  point  out  the  existence  and  usefulness  of  these 
special  distributions  and  refer  working  reliability  engineers  to  standard 
text  on  statistics-**4  for  a  more  comprehensive  discussion. 
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has  the  Student-t  distribution  with  k  degrees  of  freedom.  We  have  previously 

shown  that  p”,  i.e.,  Eq.  (30),  is  the  expectation  value  of  p  and,  therefore  the 

—  2  —2  2 

random  variable  p  -  p  has  a  zero  mean  and  variance  ou  .  Also,  since  no  fa 

2  ” 

has  been  indicated  to  have  a  x  distribution  of  (n  -  1)  degrees  of  freedom  and 
is  independently  distributed  from  p  -  p,-* 


(p  -  p)/o, 


*k  * 


M 


[ no2/ (n  -  l)o2]l/2 


(36) 


has  the  Student-t  distribution  with  (n  -  1)  degrees  of  freedom.  We  can  write 
the  confidence  level  associated  with  the  median  value  as 


(p  -  p)/o^ 


[no2/(n  -  l)o2] 1/2 


} 


(37) 


or  rearranging  terms 


_  ,  1/2 

p{p<p-t(6)n_1[^_ryI(X1-p)]  }-«  (38) 

where  the  confidence  coefficient  6  can  be  determined  from  tables  of  the 
cumulative  Student-t  distribution. 
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Returning  once  again  to  the  situation  where  we  know  o  exactly  but  not 
p;  our  best  estimate  as  to  the  failure  pdf  of  the  n  +  l  device  taken  from  the 
uniform  distribution  under  discussion  is 

f(,)  '  .-'.1/2  '’'’’I'  1  <W 

o  (2x)  u 

u 


where 


2 


2 


1 


2  2 

represents  the  summation  of  the  intrinsic  random  variation  a  and  a  . 

M 

Employing  the  Student-t  distribution,  we  can  estimate  the  failure  pdf  of  the 
n  +  1  device  by  considering 


t 


n-1 


(x  -  u)/au 

r  -2,,  2,1/2 

[na  /(n  -  l)a  ] 


(41) 


Since  the  random  variable  x  -  p  is  normally  distributed  with  zero  mean  and 
2  —2  2 

variance  oy  and  no  /o  has  a  chi-square  distribution  with  n-1  degrees  of 
freedom,  this  expression  has  the  Student-t  dependence.  We  can  estimate  the 
probability  of  the  n  +  1  device  having  a  lifetime  shorter  than  some  specified 
time  by 


P{“t(6)n-1  >  C 


(x  -  p)/o, 


tno2/(n  -  Do^] 


27172^ 


(42) 


which  on  being  rearranged  gives 

p{lnr  <  7  -  t(6)  ( 1  +  -)1/2o}  -  fi  (43) 

n+1  n—  1  n 

The  statistics  of  thermally  accelerated  aging  results  will  now  be 
analyzed.  For  this  analysis,  it  is  convenient  to  rewrite  the  failure  pdf  as 


(44) 


where  the  Arrhenius  relation  has  been  reformulated  as 

p  ■  lnt  +  AEZ 
o 


(45) 
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The  median  log  lifetime  is  a  linear  function  of  the  variable  (kT)-*,  and  the 
answer  to  our  Immediate  problem  is  provided  by  the  well-developed  simple 
linear  normal  regression  technique.  Having  provided  the  formalism  in  the  pre¬ 
ceding  .discussion,  we  may  rapidly  present  the  conclusions  of  this  analysis. 

From  Eq.  (44),  the  log-likelihood  function  is 

lni2?-  -n[lno(2ir)1/2] - -  J^-lnt  -  AEZj2  (46) 

2a 

where  the  couple  (X^.Z^)  represents  the  experimental  log  lifetime  X^,  observed 
when  the  device  is  operated  at  Z^.  The  summation  is,  of  course,  over  all 
devices  that  are  introduced  into  the  accelerated  test  program.  In  keeping 
with  our  primary  objective,  we  are  Interested  in  estimating  the  lifetime  of 
the  n  +  1  device  operated  at  a  normal  temperature  Tjj  from  the  results  of  an 
accelerated  stress  program.  Determining  the  value  of  a,  lnxo>  and  AE,  which 
maximizes  the  log  likelihood  function,  yields  the  estimators 

_  £(X.-  X)(Z  -  Z) 

AE - ± - i-: -  (47) 

I(21  -  Z)2 

lnT  -  X  -  AE  Z  (48) 

o2  -  ~  I(X±  -  1ST  -  AEZ^2  (49) 

where 


Summations  (47)  through  (51)  are  done  over  all  experimentally  observed  n 
failures.* 

Assuming  we  know  a  exactly  before  the  start  of  this  experiment,  the 
estimated  distribution  of  median  log  lifetime  is 


"  ~  .1/2  expl“  “iW  '  lnTo  ”  AEZn)  ^ 


where 


0  ,  Z2  -  2Z  Z  +  (  -  l  Z  J2 

2  2  r  N  N  '-n^i-'i 

j  -  a  [ - 5 - J 

H  l <Zj-  z)2 


The  estimated  pdf  of  the  n  +  1  device  at  ZN  is 


£<V  — 7TJ72 1-  tt(x»  -  lmo  -  SEV2] 

o  (2ir  )  2o 

u  p 


where 


2  2,2 

o  •  o  +  o„ 
u  M 


The  above  variances  can  be  much  larger  than  the  corresponding  quantity  given 
in  Eqs.  (31)  and  (40)  when  Z^  is  far  outside  the  accelerated  temperature 
range. 

An  accelerated  temperature-stress  experiment  is  not  performed  at  random 
temperatures,  but  usually  at  several  specially  selected  temperatures. 
Operating  at  several  elevated  temperatures  is  not  only  less  costly,  but  it 


For  a  normal  pdf,  the  least-squares  method  yields  Identical  results. 
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offers  the  important  advantage  of  allowing  one  to  evaluate  the  temperature- 
independence  of  the  variance.  We  will  now  consider  a  few  special  sample 
distributions.  Perhaps,  the  simplest  sample  distribution  is  a  two-temperature 
accelerated  test  with  devices  equally  distributed  between  the  two 
temperatures.  Designating  the  high-  and  low-stress  conditions  as  and  Z^, 
we  rewrite  Eqs.  (47),  (48),  (49),  (53),  and  (55)  as 


—  *L  " 

L  H 


1  f  mw  A  T9H  t  t  nO 


ln  \>  ‘  tCXL  -  *H  '  4EZL  +  4EV 


«2  - 1  n  «1  -  +  la,  -  y'l 
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2  a  „2  ^ 

°m  "  i ra  +  ° 


(56) 

(57) 

(58) 

(59) 


a 2  -  o2[l  +  -(L2  +  1)]  (60) 

u  1  n 

where  X^  and  are  the  average  log  lifetimes  at  the  low  and  high  tempera¬ 
tures.  In  these  equations,  L  represents  an  effective  experimental  lever  given 
by 


L 


E»  -  *  V 

\  -  1<ZH  +  V 


(61) 


(62) 
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and  accounts  for  the  increased  uncertainity  as  one  moves  away  from  the  accel- 

2  2 

erated  temperature  range.  Note  that  for  L  ■  o  the  variances,  o„  and  o  »are 

M  U 

equal  to  their  values  in  the  one-temperature  situation.  The  analogy  between 
this  experimental  lever  arm  and  a  similar  relation  for  the  dependence  of  the 
accuracy  of  a  rifle  on  the  ratio  of  firing  range  to  barrel  length  is  obvious. 

The  question  naturally  arises  if  an  equal  distribution  of  devices  between 
the  two  temperatures  is  statistically  the  most  effective  arrangement.  Exam¬ 
ination  of  the  above  variances  shows  them  to  be  minimized  for  the  sample 
distribution 


Uln 

h  *  TaT  +  TbT 


and,  of  course. 


n  -  n.. 


where 


A  -  1  -  L) 

B  =  1(1  +  L) 


(63) 


(64) 


(65) 

(66) 


The  minimized  variances  are 


2 


li  -  lI  +  li  +  lI)2 


(67) 


^L2 

n 


(L  >  1) 


(68) 


and 


°u  “  °2tl  L*]  (L  >  l) 


(69) 


31 


l 


From  a  practical  engineering  viewpoint,  this  minimization  is  of  little  value 
since  the  effective  lever  arm  is  typically  L  «  6. 

Another  typical  sample  distribution  is  a  three-temperature  test  in  which 
an  equal  number  of  components,  i.e.,  n/3,  are  failed  at  Z^,  ZM,  and  Z^.  Our 
estimated  parameters  are  now 


°u  "  fl2l1  +  i(T  l2+  l)J  (73) 

The  devices  failed  at  the  middle  temperature  cannot  be  used  to  estimate  the 
activation  energy.  Consequently,  for  typical  experimental  lever  arms,  the 
middle-temperature  failures  contribute  little  to  the  statistical  accuracy  of 
the  extrapolated  median  lifetime,  but  does,  of  course,  contribute  to 
estimating  o  as  indicated  by  Eq.  (49). 

A  uniform  distribution  is  of  interest  since  it  approximates  a  situation 
in  which  devices  are  stressed  at  many,  e.g.,  ten  equally  space  points  along 
the  Z  axis.  The  variance  in  this  case  is 


2-  (3L  +  1) 
n 


2 

with  a  related  expression  for  ou>  Statistically,  two-thirds  of  the  devices 
are  effectively  being  wasted  relative  to  estimating  p.  For  any  symmetric 
sample  distribution  about  Z^,  the  variance  can  be  written  as 
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(75) 


°M  “  T-^2  +  » 

where  the  D-function  is  determined  by  the  relative  sample  distribution 
throughout  the  accelerated  temperature  range.  The  D-quantity  can  be  thought 
of  as  measuring  how  effectively  the  samples  are  distributed.  The  statistical 
inefficiency,  of  course,  increases  as  the  devices  are  located  near  the  center- 
of-mass  of  the  sample.  It  is  convenient  to  define  the  statistical  efficiency 
of  an  accelerated  temperature-stress  test  as 

n  -  (DL2  +  l)"1  (76) 

The  statistical  efficiency  increases  as  both  L  and  D  are  reduced. 

The  chi-square  pdf  can  be  employed  to  estimate  the  uncertainty  in  the 
estimated  population  variance.  This  analysis  follows  the  previous  one 
presented  for  the  one-temperature  case  except  that  now 

Xn-2  *  h  2<Xi  "  IfZi)2  (77> 

a 

where  the  chi-square  distribution  now  has  n  -  2  degrees  of  freedom.  An 

additional  degree  of  freedom  is  lost,  since  at  least  two  experimental  points 

at  different  temperatures  are  needed  to  estimate  lnx  and  AE.  The  confidence 

o 

level  is  written 


(78) 


.  -  2 

In  evaluating  the  temperature-sensitivity  of  a  ,  one  should  employ  not  this 
equation,  but  Eq.  (33),  i.e.,  the  confidence  level  at  each  of  the  various 
accelerated  temperatures. 

We  will  now  repeat  the  procedure  previously  discussed  and  apply  the 
Student-t  distribution  to  determine  the  error  associated  with  »7M  and  t  . 


IP  tff" 


from  accelerated  stress  results  with  no  a  priori  knowledge  of  u  or  o.  The 
quantity 


t 


n-2 


N "  V/qM 

[no2/(n  -  2)0^] 1/2 


(79) 


has  a  Student-t  distribution,  since  -  u„  has  zero  mean  and  variance  ou, 
—2  2  N  N  n 

and  no  /o  Is  chi-square  distributed  with  (n  -  2)  degrees  of  freedom.  We, 

therefore,  can  estimate  the  probability  that  pN  will  be  lower  than  any 

particular  value  as 


(80) 


where  o  and  a  are  defined  in  Eqs.  (49)  and  (53).  Since  ou  is  proportional  to 
M  M 

o,  the  confidence  limits  for  can  be  determined.  Transforming  this 
expression  into  experimental  measured  parameters,  we  have 


<  lnTQ  +  -  t(6)n  2(o/nn)1/2}  -  5 


(81) 


where 


o2  -  — £(X  -  InT  -  4EZ.)2  (82) 

n  *  z  l  o  i 

is  our  unbiased  estimator  of  o  from  accelerated  temperature  aging. 

Our  primary  interest  is,  of  course,  in  the  probability  that  a  randomly 
selected  device  operating  at  normal  temperatures  will  have  a  lifetime  lower 
than  a  preselected  value.  This  probability  can  be  found  by  considering  the 
relation 


t 


(83) 


flntnH  ~  ^  -  AEZNJ/gu 
[no2/(n  -  2 )o]1/2 

which  can  be  shown  to  have  a  Student-t  distribution  employing  the  same  logic 
as  invoked  several  times  before.  We  write  the  desired  probability  as 

ptln,n+l  ‘  +7^]1/2)  "  «  C84) 

where  a  is  defined  in  Eq.  (82).  This  relation  expresses  the  probability  that 
t  will  be  lower  than  any  preselected  value  in  terms  of  known  experimental 
results.  It  represents  the  estimated  average  cff  that  should  be  employed  in 
estimating  the  system  loss  to  be  associated  with  this  device  and  as  such 
represents  completion  of  our  formal  statistical  analysis. 
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V.  APPLICATION  AND  SUMMARY 


To  demonstrate  the  application  of  the  preceding  analysis,  we  could  take 
thermally  accelerated  aging  results,  such  as  those  illustrated  in  Fig.  1,  and 
doggedly  estimate  the  various  probabilities  of  interest.  We  have  elected  not 
to  proceed  in  this  sterile  mode.  Instead,  we  propose  considering  a  very  arti¬ 
ficial,  but  instructive,  situation.  This  contrived  example  has  been  selected 
not  only  because  it  illustrates  the  employment  and  usefulness  of  the  preceding 
analysis,  but  also  because  it  clearly  shows  potential  pitfalls  that  one  can 
encounter  in  applying  accelerated  aging  results. 

Consider  the  problem  of  selecting  a  semiconductor  component  for  a  complex 
space  system.  Assume  that  the  system  analyst  has  determined  that  the  required 

device  should  have  an  average  failure  rate  of  less  than  10  FITs,  i.e.,  10 

failures/10^  hr,  over  the  system  useful  life  of  10^  hr  when  the  device  is 
operated  at  a  junction  temperature  of  50°C  or  less.  We  can  select  the  device 

from  six  different  lots,  all  of  which  have  been  through  an  accelerated 

temperature-stress  program.  Although  the  thermal  aging  conditions  for  each  of 
these  lots  were  different  (these  contrived  conditions  are  listed  in  Table  1), 
they  all  yielded  the  same  extrapolated  (50°C)  median  lifetime  and  logarithmic 
variance.  In  particular,  let  us  assume  that  the  extrapolated  median  lifetimes 
are  10^  hr  and  that  the  sample  distribution  logarithmic  standard 
deviations,  o  *  1.0. 

Based  on  the  preceding  analysis,  the  results  of  an  accelerated  test  pro¬ 
gram  extrapolated  to  a  given  temperature  must  be  characterized  by  four  para- 

m  A 

meters  —  not  only  and  a,  but  also  the  total  number  of  devices  in  the 
accelerated  test  and  the  statistical  efficiency  of  the  accelerated  stress  pro¬ 
gram  q.  We  have  listed  the  statistical  efficiency  of  each  experiment  in  Table 
1  with  the  parameter  o^,  which  will  prove  useful  in  discussing  this 
situation.  Both  of  these  parameters  are  dependent  on  the  device  operational 
junction  temperature,  and  the  values  listed  in  Table  2  are  for  T^  ■  50°C. 

The  probabilities  that  the  lifetime  of  the  device  during  normal  operation 
will  be  less  than  various  preselected  times,  t',  has  been  calculated  employing 
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Table  1.  Contrived  Accelerated  Temperature-Aging  Experiments 


Extrapolated 

Lot 

TH^nH 

V»L 

n 

o  /o 
u 

1 

250/16 

150/16 

5.27  x 

10~2 

1.26 

2 

250/16 

200/16 

8.64  x 

10-3 

2.15 

3 

250/8 

150/8 

5.27  x 

10-2 

1.47 

4 

250/8 

200/8 

8.64  x 

10"3 

2.87 

5 

250/4 

150/4 

5.27  x 

10-2 

1.84 

6 

250/4 

200/4 

8.64  x 

10-3 

3.93 

The  parameters 

o  /a  and  n  are 
u 

calculated 

with  application  of 

Eqs.  (60)  and  (76). 
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Table  2.  Analysis  of  Contrived  Thermal  Aging  Tests 


t'(hr) 


Lot 

6  =* 

0. 

25 

s  = 

0. 

10 

6  - 

0. 

01 

6  = 

0. 

005 

6  - 

0.0005 

0 

5.10 

X 

106 

2.77 

X 

106 

9.77 

X 

105 

7.61 

X 

105 

3.72 

X 

105 

1 

4.23 

X 

106 

1.92 

X 

106 

4.52 

X 

105 

3.13 

X 

105 

1.01 

X 

105 

2 

2.30 

X 

106 

5.98 

X 

105 

5.08 

X 

104 

2.71 

X 

104 

3.94 

X 

103 

3 

3.62 

X 

106 

1.38 

X 

106 

2.11 

X 

105 

1.25 

X 

105 

2.27 

X 

104 

4 

1.37 

X 

106 

2.11 

X 

105 

5.36 

X 

103 

1.95 

X 

103 

6.91 

X 

101 

5 

2.67 

X 

106 

7.06 

X 

105 

3.08 

X 

104 

1.09 

X 

104 

1.73 

X 

102 

6 

5.95 

X 

105 

3.49 

X 

105 

4.32 

X 

101 

4. 

71 

6.75 

X 

10~4 

The  t'  is  the  estimated  device  lifetime  at  Tj  ■  50°C  at  a  confidence  level 
(1  -  6)  (or,  equivalently,  6  failure  probabilities  at  50%  confidence  level). 
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Eq.  (84),  remembering  p  -  lnt  "  +  AEZ„  ,  and  are  listed  in  Table  2.  Also 

NON 

included  in  Table  2  (as  lot  zero)  are  the  estimated  failure  probabilities  for 
the  case  where  pN«  In  (10'  hr)  and  a  -  1.00  are  the  exact  parameters  of  a 
lognormal  failure  pdf.  For  the  reader's  convenience,  we  display  the  estimated 
average  cff  for  three  accelerated  aging  tests  as  well  as  the  cff  for  lot  zero 
in  Fig.  2. 

The  estimated  cff  of  lot  zero  can,  of  course,  be  found  directly  from  Eq. 
(16)  since  we  are  assuming  that  we  know  and  o  exactly.  One  can  calculate 
the  cff  equally  as  well  from  Eq.  (84)  by  allowing  the  sample  size  to  become 
sufficiently  large.  With  a  large  sample  size 


t(5 )  ,  +  t(6) 

n— z  ■ 


and 


(1 


+  ->) 

nn-' 


♦  1 


Under  the  condition  of  a  large  sample  size,  the  confidence  coefficient  5  ap¬ 
proaches  the  standard  normal  cff. 

7  A 

It  is  readily  apparent  that  if  p_.  -  ln(l0'  hr)  and  o  -  1.0  adequately 
characterized  the  failure  pdf,  the  failure  rate  at  times  less  than  the 
specified  useful  system  life  of  10^  hr  would  be  insignificant.  Using  the 

M  A 

sample  distributions  p„  and  o  as  the  exact  p„  and  o  is  perhaps  the  most 
serious  pitfall  to  be  associated  with  applying  accelerated  temperature-aging 
results.  For  times  less  than  p„,  the  analysis  of  the  six  contrived 
accelerated-aging  programs  indicates  an  average  probability  of  failure  much 
higher  than  one  estimates  using  p^  and  a  as  exact  failure  parameters,  l.e., 
the  lot  zero  case. 
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AVERAGE  FAILURE  PROBABILITY  (« 


TIME  (hr) 


Fig.  2.  Estimated  Average  Failure  Probability  as  Function  of  Operating 
Time  at  50°C  Inferred  From  Accelerated  Aging.  The  curve 
designated  0  is  the  cff  of  a  lognormal  pdf  with  parameters  * 
10'  hr  and  o  ■  1.00  (see  text). 


The  difference  between  the  analysis  of  lot  zero  and  Eq.  (84)  is 
fundamental  to  our  present  interest*  When  one  uses  y^  and  a  alone  to 
calculate  the  cff,  one  is  estimating  the  median  probability  of  failure  at  the 
various  operational  times  and  not  the  average  probability  of  failure  as  Eq. 

(84)  does.  At  10^  hr,  the  median  and  average  failure  probabilities  are  equal, 
but  they  diverge  as  one  goes  to  either  shorter  or  longer  times.  The 
divergence  is  a  direct  result  of  the  uncertainties  in  the  estimates 
of  y  and  o  having  a  greater  effect  at  both  short  and  long  times.  The 
probability  curves  are  symmetric  about  y^,  but  since  we  are  not  interested  in 
lifetime  greater  than  yR,  only  the  lower  portions  are  illustrated  in  Fig.  2. 

For  lots  1  and  2,  the  sample  sizes  are  sufficiently  large  that  for  6  > 
0.005,  the  error  in  the  estimated  cff  is  caused  by  the  uncertainity  in  y  and 
not  o.  In  other  words, 

t(6)3Q  -  tW)^;  6  >  0.005 

and  the  cffs  can  be  approximately  represented  by  a  lognormal  cff  with  parame- 
7  *  1/2 

ter  yN  ■  In  (10'  hr)  and  a  -  a(l  +  1/nn)  .  Although  both  lots  1  and  2  have 

failed  the  same  number  of  devices,  the  uncertainity  in  the  second  lot  at  low 
values  of  the  cff  is  much  greater  because  of  its  longer  lever  arm  or,  equiva¬ 
lently,  its  lower  statistical  efficiency. 

When  the  sample  size  is  reduced  by  one-half  (lots  3  and  4),  the  uncertain- 

—  A 

ties  in  both  y„  and  o  are  almost  equally  important  for  0.0005  <  6  <  0.01. 

N 

Above  (  ■  0.01,  the  uncertainty  in  y  dominates,  whereas  below  0.0005,  the 
uncertainty  in  o  dominates.  The  exact  region  between  where  the  uncertainty 
in  y  or  o  dominates  is,  of  course,  dependent  on  the  relative  magnitudes 
of  bW)n_2  an<*  °u/°  a8  can  clearly  seen  by  comparing  the  various  lots. 

Both  the  sample  size  and  the  statistical  efficiency  of  lots  5  and  6  are 
so  low  that  the  uncertainty  in  their  failure  rates  at  time  less  than  10^  hr 
makes  them  useless  for  our  assumed  requirement.  In  fact,  only  lot  1  should  be 


considered  as  meeting  the  requirement  of  less  than  an  average  of  10  FITs  over 
a  useful  system  lifetime  of  10^  hr. 

As  a  general  rule,  it  appears  more  appropriate  to  assume  that  the  esti¬ 
mated  average  failure  probability  follows  an  exponential  failure  pdf  with 

A 

A  “  exp(-iijj)  than  that  mn  and  a  are  the  exact  value  of  p  and  a  .  Extending 
this  suggestion  to  many  different  components,  i.e.,  components  from  different 
accelerated  aging  experiments,  should  further  Increase  the  accuracy  of  this 
approximation.  Considering  that  it  is  not  difficult  to  calculate  an  effective 
failure  rate  with  the  foregoing  analysis  when  the  accelerated  aging  conditions 
are  adequately  characterized,  the  assumption  of  an  exponential  pdf  should  be 
considered  only  when  the  thermal  stress  results  are  not  properly  characterized. 

The  results  of  this  analysis  are  best  summarized  by  emphasizing  that  any 
accelerated  temperature-aging  test  should  be  characterized  by  more  parameters 
than  just  the  logarithmic  variance  of  the  failed  sample  distribution  and  the 
extrapolated  median  lifetime  (or,  since  all  users  will  not  operate  at  the  same 
junction  temperature,  Tq  and  Af).  This  is  especially  important  when  one 
wishes  to  estimate  the  failure  rate  at  times  much  different  than  the  extrapolated 
median  lifetime.  To  properly  estimate  the  average  probability  of  failure  at 
times  much  less  than  pN>  the  total  number  of  devices  in  the  accelerated 
temperature-aging  program  should  be  provided  as  well  as  each  failure  tempera¬ 
ture.  These  parameters  are  then  used  to  calculate  the  probability  of  failures 
at  low-level  cff,  and  an  estimated  average  failure  rate  is  calculated.  If  the 
average  failure  cff  meets  the  system  requirements,  one  is  justified  in  con¬ 
cluding  that  the  statistical  requirement  has  been  fulfilled  and  one  should 
then  turn  to  the  consideration  of  potential  systematic  errors  in  applying 
accelerated  temperature-aging  results  —  a*  subject  beyond  our  limited 
objective. 
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ABBREVIATIONS  AND  SYMBOLS 

Cumulative  failure  function 

Cumulative  failure  function  (cff)  of  random  variable  (x) 

Failure  probability  density  function  (pdf)  of  random  variable  (x) 
Boltzmann’s  constant 

(ZN  ”  ZM^ZL  “  hO 

Natural  logarithm 
Logarithm  to  base  10 

Total  sample  size  in  accelerated  temperature  stress  (all  failed) 
Sample  size  at  T^  (all  failed) 

Probability  density  function 

Absolute  "junction"  temperature 

Highest  accelerated  temperature 

Lowest  accelerated  temperature 

Normal,  i.e. ,  desired,  operational  temperature 

Time,  random  variable 

Student-t  distribution  with  k  degrees  of  freedom  associated  with 
confidence  coefficient  6 

In 

In  t,  random  variable 

(kV1 

(ZH  +  Zl)/2 

Activation  energy  of  Arrhenius-dependent  failure  mechanism 
Confidence  coefficient 
Device  failure  rate 


O  (VI  M 


u 

“i 

**N 

2 

a 

-2 

a 

“2 

o 


T 


Median  llfetlme/log  lifetime 

Mean  of  sample  distribution  at 

Mean  of  sample  distribution  extrapolated  to 

Variance/log  variance 
2 

Estimator  of  o 

2 

Unbiased  estimator  of  a 
'’Experimental"  failure  time  of  device 

Pre-exponential  factor  of  Arrhenius-dependent  failure  mechanism 
Chi-squares  distribution  with  k  degrees  of  freedom 
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LABORATORY  OPERATIONS 


The  Laboratory  Operations  of  The  Aerospace  Corporation  le  conducting 
experimental  and  theoretical  Investigations  naceeeary  for  Che  evaluation  and 
application  of  eclentlflc  advancaa  to  new  military  concepts  and  eystame.  Var- 
aatlllty  and  flexibility  have  been  developed  to  a  high  degree  by  the  laborato¬ 
ry  personnel  In  dealing  with  the  nany  problems  encountered  In  the  Nation's 
rapidly  developing  space  ay a tana.  Expertlae  In  the  latest  scientific  develop¬ 
ments  Is  vital  to  the  accomplishment  of  tasks  related  to  these  problems.  The 
laboratories  that  contribute  to  this  research  are: 

Aarophyslca  Laboratory:  Aerodynamics;  fluid  dynamics;  plaamadynamlcs; 
chemical  kinetics;  engineering  nachaalee;  flight  dynamics;  heat  transfer; 
high- power  gas  lasers,  continuous  and  pulsed,  11,  visible,  BV;  laser  physics; 
laser  resonator  optica;  laser  effects  and  countermeasures . 

Chemistry  and  Physics  Laboratory:  Atmospheric  reactions  and  optical  back- 
grounds;  radiative  transfer  and  atmospheric  transmission;  thermal  and  state- 
specific  reaction  rates  In  rocket  plumes;  chemical  thermodynamics  and  propul¬ 
sion  chemistry;  laser  Isotope  separation;  chemistry  aad  physics  of  particles; 
apace  environmental  and  contamination  affects  aa  spacecraft  materials;  lubrica¬ 
tion;  surface  chemistry  of  Insulators  ami  conductors;  cathede  materials;  sen¬ 
sor  materials  aad  sensor  optics;  spplled  laser  spectroscopy;  atomic  frequency 
standards i  pollution  aad  toxic  materials  monitoring. 

Electronics  Research  Laboratory;  Electromagnetic  theory  aad  propagation 
phenomena;  microwave  aad  semiconductor  devices  aad  latagratad  circuits;  gusa- 
tun  electronics,  lasers,  and  electro-optics;  communication  eel sac as,  applied 
electronics,  superconducting  and  electronic  device  physical  allllaeter-anva 
and  far-lafrsrad  technology. 

Materials  Sciences  Laboratory;  Development  of  new  eater la la;  composite 
materials;  graphite  and  ceramics;  polymeric  materials;  weapons  affects  aad 
hardened  materials;  materials  for  electronic  devices;  dimensionally  stable 
materials;  chemical  aad  structural  analyses;  stress  corrosion;  fatigue  of 
aatala. 

Specs  Sciences  Laboratory:  Atmospheric  end  ionospheric  physics,  radia¬ 
tion  from  the  ataoaphsra,  density  and  composition  of  the  atmosphere,  auroras 
and  alrglow;  magnetos pharlc  physics,  cosmic  rays,  generation  and  propagation 
of  plasma  waves  in  the  aagnstos phare;  solar  physics,  x-ray  astronomy;  the  effect! 
of  mtclaar  explosions,  aagaetlc  stoma,  and  solar  activity  on  the  earth's 
atmosphere,  ionosphere,  and  nagnatoaphara;  the  effects  of  optical,  electronag-  . 
ns tic,  aad  particulate  radiations  in  apnea  on  apnea  systems. 


